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Results are reported on the measurement of the atmospheric neutrino-induced muon flux at 
a depth of 2 kilometers below the Earth's surface from 1229 days of operation of the Sudbury 
Neutrino Observatory (SNO). By measuring the flux of through-going muons as a function of 
zenith angle, the SNO experiment can distinguish between the oscillated and un-oscillated portion 
of the neutrino flux. A total of 514 muon-like events are measured between —1 <cosSzonith< 0.4 in 
a total exposure of 2.30 x 10^" cm^ s. The measured flux normalization is 1.22 ± 0.09 times the 
Bartol three-dimensional flux prediction. This is the flrst measurement of the neutrino-induced 
flux where neutrino osciUations are minimized. The zenith distribution is consistent with previ- 
ously measured atmospheric neutrino oscillation parameters. The cosmic ray muon flux at SNO 
with zenith angle cos6'zonith> 0.4 is measured to be (3.31±0.01 (stat.)±0.09 (sys.)) x 10~^" /i/s/cm^. 



PACS numbers: 14.60.Lm, 96.50.S-, 14.60.Pq 
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I. INTRODUCTION 



II. THE SUDBURY NEUTRINO 
OBSERVATORY 



Atmospheric neutrinos are produced from the decay of 
charged mesons created by the interactions of primary 
cosmic rays with the Earth's atmosphere. Atmospheric 
neutrinos can be detected either via direct interactions 
within the fiducial volume of a given detector or indi- 
rectly from the observation of high-energy muons cre- 
ated via the charged current interaction 1^^^ + N + X 
on materials that surround the detector. Although the 
latter process produces muons propagating at all zenith 
angles, overhead portions of the sky are typically domi- 
nated by cosmic-ray muons created in the Earth's upper 
atmosphere. 



The flux of atmospheric neutrinos has been a topic of 
study since the mid-1960's. Early experiments [H [2] in- 
ferred the presence of atmospheric neutrinos by measur- 
ing the muon flux created by neutrino interactions taking 
place in rock surrounding a detector. Subsequent studies 
of the atmospheric neutrino flux as a function of zenith 
angle [3, 4, 5], the ratio of electron and muon neutri- 
nos 1^ lL i8j i9j [ID] , and combined measurements have 
provided a more direct measurement of the atmospheric 
neutrino flux and revealed evidence for neutrino oscil- 
lations. Results gathered from these experiments have 
been further verified by long baseline accelerator mea- 
surements |12[ 113] , thereby providing strong constraints 
on the neutrino oscillation parameters. 

The Sudbury Neutrino Observatory (SNO) is located 
in the Vale-Inco Creighton mine in Ontario, Canada at 
a depth of 2.092 km (5890 ± 94 meters water equiva- 
lent) with a flat overburden [14'. The combination of large 
depth and flat overburden attenuates almost all cosmic- 
ray muons entering the detector at zenith angles less than 
cos0zcnith= 0.4. Bccausc of this depth, SNO is sensitive 
to neutrino-induced through-going muons over a large 
range of zenith angles, including angles above the hori- 
zon. 



This paper presents a measurement of the flux of 
muons traversing the SNO detector. Measuring the 
through-going muon flux, as a function of zenith angle, 
for cos 6'zonith< 0.4 provides sensitivity to both the oscil- 
lated and un-oscillated portions of the atmospheric neu- 
trino flux. Measuring the muon angular spectrum above 
this cutoff provides access to the flux of cosmic-ray muons 
created in the upper atmosphere. This paper is divided as 
follows: Section |ll] describes the experimental details of 
the SNO detector, Section |ITT] describes the Monte Carlo 
model used to predict the observed muon flux. Section IV 
describes the data collection and event reconstruction, 
and Section |V] discusses the signal extraction and error 
analysis used for the measurements presented herein. 



The Sudbury Neutrino Observatory is located 
at 46°28'30"N latitude (56° 33' magnetic north), 
81°12'04"W longitude near the city of Sudbury, Ontario. 
The center of the detector is at a depth of 2092 ± 6 me- 
ters from the Earth's surface. The Earth's surface im- 
mediately above the SNO detector is 309 meters above 
sea level. Within a 5 km radius on the surface above 
the detector, the local topology lies between 300 and 320 
meters above sea level, although small localizations with 
±50 meters level variations do occur. The norite rock 
that dominates the overburden is mostly oxygen (45%), 
silicon (26%), aluminum (9%), and iron (4%). A com- 
bination of bore samples taken at different depths and 
gravity measurements taken at the surface show varia- 
tions in the rock density, from 2.8 g/cm^ near the surface 
and closer to 2.9 g/cm'^ in the vicinity of the detector. 
A fault line located 70 meters southwest of the detector 
serves as a boundary to a deposit of granite/gabbro rock 
of similar density (2.83 ± 0.10 g/cm'^) but slightly differ- 
ent chemical composition {{Z'^/A) of 5.84 versus 6.01). 
An average rock density of 2.83 ± 0.05 g/cm'^ is used in 
overburden calculations, independent of depth. The un- 
certainty in the density takes into account the variation 
as measured in the rock volume surrounding the detector. 
The total depth to the center of the SNO detector, tak- 
ing into account air and water flUed cavities, is 5890 ±94 
meters water equivalent. 

The SNO detector itself includes a 600.5 cm radius 
acrylic vessel flUcd with 99.92% isotopically pure heavy 
water (D2O). The 5.5-cm thick acrylic vessel is sur- 
rounded by 7.4 kilotons of ultra-pure H2O encased within 
an approximately barrel-shaped cavity measuring 34 m in 
height and 22 m (maximum) in diameter. A 17.8-meter 
diameter stainless steel geodesic structure surrounds the 
acrylic vessel. The geodesic is equipped with 9456 20- 
cm photo-multiplier tubes (PMTs) pointed toward the 
center of the detector. A non-imaging light concentrator 
is mounted on each PMT to increase the total effective 
photocathode coverage to 54%. 

SNO is primarily designed to measure the solar neu- 
trino flux originating from ®B decay in the sun above 
a threshold of several MeV by comparing the observed 
rates of the following three reactions: 



i^x + er + e'^ (ES) 

Ve + d e- +P + P (CC) 
Vx + d er +p + n (NC) 



(1) 



The charged current (CC), neutral current (NC), and 
elastic scattering (ES) reactions outlined above are sensi- 
tive to different neutrino flavors. Data taking in the SNO 
experiment is subdivided into three distinct phases, with 
each phase providing a unique tag for the final states of 
the neutral current interaction. In the first phase, the 
experiment ran with pure D2O only. The neutral current 
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reaction was observed by detecting the 6.25-MeV 7-ray 
following the capture of the neutron by a deuteron. For 
the second phase of data taking, approximately 0.2% by 
weight of purified NaCl was added to the D2O to enhance 
the sensitivity to neutrons via their capture on "^^Cl. In 
the third and final phase of the experiment, 40 discrete 
■^He or *He-filled proportional tubes were inserted within 
the fiducial volume of the detector to enhance the capture 
cross-section and make an independent measurement of 
neutrons by observing their capture on ^He in the pro- 
portional counters. Results from the measurements of the 
solar neutrino flux for these phases have been reported 
elsewhere [13 [H HZl UHl US] . 

Muons entering the detector produce Cherenkov light 
at an angle of 42° with respect to the direction of the 
muon track. Cherenkov light and light from delta rays 
produced collinear to the muon track illuminate an aver- 
age of 5500 PMTs. The charge and timing distribution 
of the PMTs is recorded. The amplitude and timing re- 
sponse of the PMTs is calibrated in situ using a light dif- 
fusing sphere illuminated by a laser at six distinct wave- 
lengths [20] . This "laser ball" calibration is of particular 
relevance to the muon analysis since it provides a tim- 
ing and charge calibration for the PMTs which accounts 
for multiple photon strikes on a single PMT. Other cal- 
ibration sources used in SNO are described in the refer- 
ences [T1I5T]. 

For a period at the end of the third phase of the exper- 
iment, a series of instrumented wire tracking chambers 
and scintillator panels were installed immediately above 
the SNO water cavity to provide a cross-check on the ac- 
curacy of the muon reconstruction algorithm. Details of 
the apparatus and results obtained from this calibration 
are reported later in this paper. 



III. SIMULATION OF SIGNAL AND 
BACKGROUND EVENTS 

Candidate neutrino-induced through-going muon 
events can arise from a variety of sources. These include: 
(a) muons created from neutrino-induced interactions 
in the rock surrounding the SNO cavity; (b) muons 
created from neutrino-induced interactions in the H2O 
volume surrounding the PMT support structure; (c) z^^ 
interactions that take place inside the fiducial volume 
but are misidentified as through-going muons; (d) i^e 
interactions that take place inside the outlined fiducial 
volume but mis-reconstruct as through-going muons; (e) 
cosmic-ray muons created in the upper atmosphere that 
pass the zenith angle cut; and (f) events created by in- 
strumental activity in the detector. The first three event 
types are proportional to the i'^ atmospheric neutrino 
flux and can undergo oscillations. The i^e-induced flux 
is also proportional to the overall atmospheric neutrino 
flux, but the currently measured neutrino oscillation 
parameters indicate that their probability for undergoing 
oscillations is highly suppressed. The last two entries 



constitute a genuine source of background to the signal. 

In order to understand the measured neutrino-induced 
flux, a proper model of the initial neutrino flux and sub- 
sequent propagation is necessary. SNO uses the Bartol 
group's three-dimensional calculation of the atmospheric 
neutrino flux [22]. Figure [l] shows the predicted flux 
for cosmic rays and muons from the interaction of muon 
neutrinos and anti-neutrinos as a function of muon en- 
ergy. The neutrino energy spectrum is correlated with 
the primary H and He cosmic flux, both of which are 
strongly constrained by data. The uncertainties that 
dominate the neutrino energy distribution relate to the 
primary cosmic-ray energy spectrum, the tt and K pro- 
duction ratio, and hadronic cross sections. Treatment of 
the systematic errors in the neutrino flux is discussed in 
greater detail in Ref [25] and [Si]. Current estimates of 
the neutrino flux uncertainties are approximately ±15% 
and depend strongly on neutrino energy. Because the 
normalization of the neutrino flux and the energy spec- 
tral shape are highly correlated, the flts to the data re- 
ported herein assume a fixed neutrino energy spectrum. 
We also assume that the flux and energy spectra do not 
change significantly with solar activity. Although vari- 
ations throughout the solar cycle are expected, the ma- 
jority of this variation is confined to neutrinos of energy 
below 10 GeV, so the impact on the fluxes predicted at 
SNO is expected to be small. A flux uncertainty of ±1% 
is included to account for variations due to solar cycle 
activity. 




4.5 



FIG. 1: The muon energy spectrum (given in logio(-E'/j)) in the 
SNO detector from cosmic-ray muons (triangles) as predicted 
from MUSIC, and from neutrino-induced muons (boxes) cre- 
ated in the surrounding rock as predicted from the Bartol 
3D Monte Carlo. The expected spectrum after oscillations is 
shown by the solid line. The distributions are not to scale. 

Neutrino interactions in the rock surrounding the de- 
tector are simulated by the NUANCE v3 Monte Carlo 
neutrino event generator [25] . NUANCE includes a com- 
prehensive model of neutrino cross sections applicable 
across a wide range of neutrino energies. Neutrino quasi- 
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elastic interactions are modeled according to the formal- 
ism of Llewellyn-Smith |26]. A relativistic Fermi gas 
model by Smith and Moniz [27] is used to model the low 
momentum transfer effects in the nucleus. The quasi- 
elastic cross-section depends strongly on the value of the 
axial mass used in the axial form factor. Recent mea- 
surements from K2K and MiniBooNE show a 
higher value of the axial mass than previously reported 
(maxiai = 1-20 ± 0.12 GeV and niaxiai = 1-23 ± 0.20, re- 
spectively) . Though this analysis uses the previous world 
average for the axial mass (maxiai — 1.03±0.15 GeV) [3D], 
the systematic uncertainty encompasses these more re- 
cent measurements. For the Fermi gas model, we assume 
a Fermi momentum of 225 MeV/c and a binding energy 
of 27 MeV for light elements such as oxygen, carbon, and 
silicon, and a 2.22 MeV binding energy for deuterium. 

For neutrino reactions where a single charged or neu- 
tral pion is resonantly produced, NUANCE employs a 
modified model of Rein and Sehgal |31| . Experimental 
constraints on this cross-section are of order ±20%, not as 
strong as those placed on the quasi-elastic cross-section. 
Many of the other parameters used for the quasi-elastic 
cross-section are also used for this process. 

The largest contributor to the atmospheric neutrino- 
induced muon flux is deep inelastic scattering of neutri- 
nos in the surrounding rock where the hadronic invariant 
mass is above 2 GcV/c^. The uncertainty on the cross- 
section for this process is strongly constrained by ac- 
celerator produced high energy neutrino experiments to 
±3% [321 221 125 ■ The transition between resonance and 
deep inelastic scattering cross-section uses the method- 
ology developed by Yang and Bodek [3S]. Other minor 
processes that can produce muons in the final states, such 
as coherent pion production and — e~ scattering, are 
also included. 

Transport of muons through the rock from neutrino- 
induced interactions is calculated using the PROPMU 



muon transport code f36|, which is integrated into the 
NUANCE Monte Carlo framework. Rock compositions 
and densities consistent with measured values are imple- 
mented in PROPMU. Simulation of muon transport in 
the D2O and H2O and subsequent detector response is 
handled by the SNO Monte Carlo and Analysis (SNO- 
MAN) code. SNOMAN propagates the primary par- 
ticles and any secondary particles (such as Compton 
electrons) that are created, models the detection of the 
optical photons by the PMTs, and simulates the elec- 
tronics response. The SNOMAN code has been bench- 
marked against calibration neutron, gamma, and electron 
data taken during the lifetime of SNO. With the excep- 
tion of a few physics processes (such as optical photon 
propagation), widely used packages such as EGS4 [37], 
MCNP [38^ and FLUKA [3S] are used in SNOMAN. Ex- 
plicit muon energy loss mechanisms such as ionization, 
pair production, bremsstrahlung, muon capture and de- 
cay, and photonuclear interactions are all included in the 
simulation, allowing modeling of the muon track from 
rest energies up to several TeV. Energy losses due to 
photonuclear interactions are simulated using the for- 
malism of Bezzukov and Bugaev [IDIHI]- Production of 
secondary particles from muon interactions, which con- 
tributes to the total energy deposited in the detector, is 
included in the model as described above. 

In addition to the through-going signal from atmo- 
spheric neutrinos, a number of backgrounds which have 
muon signatures are simulated in the analysis. These in- 
clude t'^ interactions inside the H2O and D2O volumes 
of the detector, and Ve interactions that either have a 
muon in the final state or are misidentified as a through- 
going muon. Cosmic-ray muons incident on the detector 
constitute an additional source of high-energy muons and 
are treated separately. Their flux is estimated using the 
formalism of Gaisser [42], 



-^^\^Io-E-^ ^4 ^ + , ) cm-2 sr-i GeV-i (2) 

^ ^ 115 GcV ^ 850 GcV 



where and 0^ are the muon energy and zenith an- 
gle at the Earth's surface, 7 = 2.77 ± 0.03 is the muon 
spectral index, and Iq is a normalization constant. Al- 
though Eq. [2]is inaccurate at low energies, the minimum 
energy required for surface muons to reach the SNO de- 
tector is ~ 3 TeV. Transport of such high-energy muons 
in the rock is performed by the MUSIC muon transport 
code [15] . The average energy of these muons as they en- 
ter the SNO detector is ~ 350 GeV. After incorporating 
the detector response, simulated events for cosmic-ray 
and neutrino-induced muon candidates are used to con- 
struct probability distribution functions (PDFs). The 



reconstructed zenith angle is used to establish PDFs for 
both signal and background. 



IV. EVENT SELECTION AND 
RECONSTRUCTION 

A. Livetime 

The data included in this analysis were collected dur- 
ing all three SNO operation phases. During the initial 
phase, data were collected from November 2, 1999 until 
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May 28, 2001. During the second phase of SNO, data 
were recorded between July 26, 2001 and August 28, 
2003. The third and final phase of SNO operations col- 
lected data between November 27, 2004 and November 
28, 2006. Data were collected in discrete time intervals, 
or runs, that range from 30 minutes to 96 hours in length. 
Runs that were flagged with unusual circumstances (pres- 
ence of a calibration source in the detector, maintenance, 
etc.) were removed from the analysis. The raw live time 
of the data set is calculated using a GPS-sychronized 10 
MHz clock on a run- by-run basis, checked against an in- 
dependent 50 MHz system clock, and corrected for time 
removed by certain data selection cuts. The livetime of 
the dataset here is 1229.30±0.03 days. The livetime used 
in this analysis differs from previously published analy- 
ses because requirements that have a strong impact on 
solar neutrino analyses, such as radon activity levels, are 
relaxed here. 



B. Event Reconstruction 

The SNO detector has a nearly ideal spherically sym- 
metric fiducial volume. The algorithm used to recon- 
struct muon candidates makes use of this symmetry in 
finding the best fit track. A two-tiered algorithm is used 
whereby a preliminary track is reconstructed which later 
serves as a seed for a more comprehensive fit to the 
muon candidate event. In the preliminary fit, the en- 
trance position is determined by looking at the earliest 
hit PMTs, and the exit position is determined by the 
charge- weighted position of all fired PMTs. In our spher- 
ical geometry, the impact parameter is the distance from 
the center of the sphere to the midpoint of the line con- 
necting the entrance and exit points. The fitter corrects 
the track fit for biases in charge collection and geome- 
try, and provides a first estimate of the direction of the 
incoming muon. The first-order reconstructed track is 
then passed to a full likelihood fit to determine the muon 
track parameters to greater accuracy. The likelihood fit 
uses three distributions: (a) the number of detected pho- 
toelectrons, (b) the PMT charge distribution, and (c) the 
PMT timing distribution. The charge and timing distri- 
butions are conditional on the number of photoelectrons 
incident on a given PMT. These distributions are cor- 
rected according to biases measured during laser calibra- 
tions. Figure [2] shows the timing distribution expected 
for different numbers of photoelectrons that are above de- 
tection threshold. Both changes in pre and post-pulsing 
and biasing can be seen in the timing distributions due 
to multiple photon hits on a given PMT event. 

The use of conditional distributions helps remove re- 
construction biases due to multiple photoelectrons de- 
tected on a single PMT. This is important for impact 
parameter values close to the PMT support structure 
(6p ~ 830 cm) . The quality of reconstruction for the bias 
of the fitted track and for the mis-reconstruction angle as 
a function of impact parameter were examined. Figure [3] 




Time (ns) 

FIG. 2: The probability distribution of PMT firing times 
based on simulation for events with one (dashed), two (dot- 
ted), and three (solid) photons striking the photocathode. 



shows the cosine of the mis-reconstruction angle, defined 
as the dot-product between the true muon direction vec- 
tor and the reconstructed vector. Approximately 87% 
(97%) of all simulated muons with an impact parameter 
of less than 830 cm reconstruct within 1° (2°) of the true 
track direction; respectively (see Figure |3|. Monte Carlo 
studies also show that bias effects on the reconstructed 
impact parameter to be less than ±4 cm (see Figure |4]). 

C. Event Selection 

After run selection, low-level cuts are applied to mea- 
surements of PMT outputs before reconstruction in order 
to separate through-going muon candidate events from 
instrumental background activity. We require a mini- 
mum of 500 valid (or calibrated) PMT hits for an event 
to be a muon candidate. Events with more than 250 
hit PMTs within a 5 /is window of a previously tagged 
event, or when 4 or more such events occur within a 2 s 
window, are identified as burst events. Burst events are 
often associated with instrumental backgrounds and are 
removed from the analysis. Instrumental activity typi- 
cally has broad PMT timing distributions and/or low to- 
tal charge; events with these characteristics are removed. 
Finally, events that possess 4 or more hit PMTs in the 
aperture of the D2O vessel (neck) are removed from the 
data to eliminate occurrences where light enters the de- 
tector from the top of the acrylic vessel. 

A series of high-level analysis cuts use reconstructed 
track parameters to isolate a pure through-going muon 
data set. A cut on the reconstructed impact parameter of 
6p < 830 cm is applied to the data to ensure accurate re- 
construction of through-going tracks. These cuts define 
a total fiducial area of 216.42 m^ and a minimal track 
length of 367 cm. The minimum (mean) muon energy 
needed to traverse this length of track is 0.8 (2.6) GeV. 
Muon events characteristically produce ample amounts 
of light in the detector. The number of Cherenkov pho- 



0.5 



1.5 



2.5 



Q I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



mis-reconstrucl 



(degrees) 



FIG. 3: Distribution of mis-reconstructed zenith angle (left) and as a function of normalized impact parameter, = 
6^/(850 cm)^, (right) for Monte Carlo cosmic-ray muons. The points and error bars on the right-hand side plot refer to 
the mean and RMS of the mis-reconstructed angle, respectively. 



tons produced by the muon, scaled by the appropriate 
detection efficiency for photons produced at the given 
impact parameter, is reconstructed for each candidate 
event. Each track is required to possess a minimum of 
2000 detected photoelectrons. A cut is also made on the 
estimated energy loss (dE/dX) of the muon. The quan- 
tity dE/dX is determined from the amount of detected 
light, corrected for geometric and photon attenuation ef- 
fects, divided by the reconstructed track length. The 
dE/dX variable depends on both the ionization and ra- 
diation losses, and has a peak at around 225 MeV/m. 
Events with dE/dX > 200 MeV/m are retained for fur- 
ther analysis. Further, a cut is imposed on the fraction 
of photoelectrons within the predicted Cherenkov cone 
for the muon track, and on the timing of these in-cone 
photons. Finally, a linear combination (Fisher discrim- 
inant) formed from the fraction of in-time hits and the 
time residuals from the muon fit is used to reduce the 
contamination of contained atmospheric neutrino events 
in our final data sample. A list of all cuts and their effects 
on the data is shown in Table |T] and in Figure |6] 

The reconstructed cosmic-ray tracks, after all selec- 



tion criteria are applied, exhibit a flat distribution ver- 
sus impact area, as expected (see Figures |4]and|5]). The 
reconstruction efficiency is also robust to a number of 
changes in the optical and energy loss model of the recon- 
struction. Monte Carlo simulation shows that changes in 
the parameters of the detector model, including Rayleigh 
scattering, secondary electron production, PMT photo- 
cathode efficiency, and PMT angular response, all have 
minimal impact on the reconstruction performance. An 
uncertainty on the reconstruction efficiency of ±0.3% is 
assigned due to detector model dependence. The effi- 
ciency of the event selection depends most sensitively 
on the energy loss parameter, dE/dX. We determine 
the energy loss model uncertainty on the reconstruc- 
tion efficiency by studying the level of data-Monte Carlo 
agreement. For through-going cosmic-ray muons it is 
±0.2% and for neutrino-induced muons it is ±2.5%. Sim- 
ilar uncertainties in reconstruction efficiency arise from 
the PMT charge model invoked in reconstructing events 
and in the rejection of events from the linear discrimi- 
nat cut previously mentioned. This leads to a ±0.05% 
(±1.0%) and ±0.37% (±2.1%) uncertainty on the cosmic- 
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FIG. 4: Left: Comparison of true and reconstructed impact parameter versus normalized impact parameter, — &^/(850 cm)^ 
for Monte Carlo cosmic-ray muons. Data points indicate mean and error bars for a given impact parameter. Right: Projection 
of difference between reconstructed and generated tracks. Dashed lines indicate uncertainty in impact parameter reconstruction 
as adopted for this analysis. 
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FIG. 5: Distribution of the normalized impact parameter, 
= ^'^/(850 cm)^ for cosmic-ray muon data (points) and 
Monte Carlo (solid line). 



ray (neutrino-induced) flux from the charge and hnear 
discriminant cuts, respectively. The differences seen in 
the two muon sources are due to the differences in muon 
energy distribution. Monte Carlo studies of cosmic-ray 
events that pass through the SNO detector show the to- 
tal event selection cut efficiency to be 99.2% for through- 
going muons. 



D. Quality Checks and Calibration 

Neutrino-induced muons have a minimum energy of 
about 2 GeV with significant intensity extending into 
the hundreds of GeV range. There is no readily available 
controlled calibration source that can provide multi-GeV 
muons as a benchmark to test the reconstruction algo- 
rithms. Instead, a number of checks have been carried 
out to test the performance of the Monte Carlo by com- 
paring with data. 

The majority of checks are performed using muons that 
reconstruct in the downward direction (cos (^zenith) > 
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TABLE I; Summary of low- and high-level cleaning cuts ap- 
plied to the data and their effect on the data population. Cuts 
are applied in sequence as they appear in this table. 



Level 


Type of Cut 


No. of Events 




Raw number of tubes hring > 250 


378219 




Timing and burst requirements 


375374 


Low 


Number of calibrated tubes firing 


100396 




Raw PMT charge requirement 


85703 




Raw PMT tmimg RMS 


84414 




TVT 1 r 111 r* • 

JNumber of neck tubes farmg 


84038 




Impact parameter < 830 cm 


80165 




Fit number of photoelectrons 


79998 


High 


Energy loss (dE/dX) 


79268 




Linear discriminant cut 


77321 




Cherenkov cone in-time fraction 


77321 




Cherenkov cone fraction of tubes firing 


77263 


Zenith 


COs6lzemth> 0.4 


76749 




COs6'zenith< 0.4 


514 




2.4 2.6 2.8 3 3.2 3.4 3.6 3.8 4 

Log (Nhits) 



FIG. 6: The number of muon candidate events as a function 
of the log of the number of tubes that fire. Plot shows events 
with no cuts applied (solid), after the burst cut (dashed), cal- 
ibrated tubes cut (dotted), all low-level cuts (red) and high- 
level cuts (filled area) are applied. 

0.4). Although the total energy of these cosmic-ray 
muons extends much higher than those from neutrino- 
induced muons, the amount of energy deposited in the 
detector by both is very similar. The slight differences 



between the two energy loss distributions arise mainly 
from the relative fraction of radiative processes contribut- 
ing to the energy loss mechanisms. Figure [7] shows the 
distribution of total hit PMTs (NHit) for reconstructed 
cosmic-ray muons. In general, there is good agreement 
between data and Monte Carlo simulations. 
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FIG. 7: The distribution of hit tubes (NHit) for muon events 
that pass all cuts for data (crosses) and Monte Carlo (line). 
The Monte Carlo has been normalized to the total number of 
events seen in the data. Only statistical errors are shown in 
the figure. 

The neutrino oscillation analysis is particularly sen- 
sitive to two parameters: (a) the fitter bias in recon- 
structing events at the edge of the impact parameter 
acceptance and (b) the angular resolution of the muon 
zenith angle. The former affects the fiducial area of the 
experiment, while the latter affects the neutrino angu- 
lar distribution, thereby affecting the oscillation param- 
eter extraction. To test the accuracy of the muon track 
reconstruction, data and Monte Carlo distributions for 
cosmic-ray muons are compared at high impact parame- 
ters. A chi-square (Axi) test is performed between data 
and Monte Carlo simulations for impact parameter dis- 
tributions under different models of impact parameter 
bias and resolution: (a) a constant shift (6^ — bfj_ + 5x), 
(b), a linear bias {b'^ = 6^ • (1 -f Sx)), and (c) a larger 
impact parameter resolution. Results are summarized 
in Table |TT] and show that our reconstruction model is 
consistent with the data and a small reconstruction bias 
(% < 4 cm at 68% C.L.). 

The muon fitter uses both time and charge to recon- 
struct muon tracks. To test the robustness of the algo- 
rithm, tracks are fit under two conditions, using "charge 
only" and "time-only" information in order to search for 
potential biases in reconstruction (Figure |8| . Differences 
between the two tracking methods are well-modeled by 
the Monte Carlo simulations (see also Table . The ob- 
served 1 . 1 cm shift is interpreted as the lower limit on the 
accuracy of the muon impact parameter reconstruction. 

A more comprehensive test of the muon tracking al- 
gorithm is to compare tracks reconstructed in the SNO 
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FIG. 8: Dilference in reconstructed impact parameter (b^) 
using "time-only" and "charge-only" information in the like- 
lihood minimization scheme for muons that pass all analysis 
cuts (points). The mean difference between the two meth- 
ods shows a 1.1 cm offset in the reconstruction of the impact 
parameter in comparison to simulations (solid line). 



detector with an external charged particle tracking sys- 
tem. A muon tracker was installed immediately above 
the SNO detector. The apparatus took data for a period 
at the end of the third phase of the SNO experiment. A 
total of four wire planes, each spanning an area of ap- 
proximately 2.5 X 2.4 m^ and containing 32 instrumented 
wire cells, was arranged in alternating orthogonal coor- 
dinates to provide two dimensional track reconstruction. 
The overall structure was tilted at a 54-degree angle due 
to pre-existing space constraints. Three large scintillator 
panels spanned the area covered by the wire chambers, 
and provided the trigger for the wire chamber readout. 
A common trigger was also sent to the main SNO data 
acquisition system to synchronize events. A total of 94.6 
days of live time was recorded by the muon calibration 
unit. Track reconstruction of muon candidates from this 
instrument were compared with the SNO tracking algo- 
rithm. Further details on this calibration system will be 
described in a future article. 

High-energy muons passing through this external 
muon tracking system and the SNO detector were recon- 
structed by both systems, providing a calibration check 
for the main SNO tracking algorithm. A total of 30 
tracks were used for comparisons of track reconstruc- 
tion between both systems. This test confirms the ac- 
curacy of the muon track reconstruction to better than 
0.62° ± 0.12°. A small shift in the reconstructed impact 
parameter is observed in the data, consistent with the 
limits from the previously mentioned tests. We take the 



uncertainty on the impact parameter reconstruction as 
±4 cm. 

TABLE IL List of consistency checks for the through-going 
muon analysis, including AxS tests on the radial distribution 
of cosmic-ray muons, charge and time reconstruction differ- 
ences, and the external muon tracking data. See text for more 
details. 



Method 


Bias 


Axl Tests 


Impact Parameter Bias 


+ 1.0% 
-0.6% 


Impact Parameter Shift 


tli cm 


Impact Parameter Resolution 


+8.5 cm 


Charge-Time Reconstruction 


Impact Parameter Shift 


±1.1 cm 


External Muon Chambers 


Angular Resolution 


< 0.6° 


Impact Parameter Bias 


4.2 ±3.7 cm 



A final check is performed on the time interval dis- 
tribution between muon events. A fit to an exponential 
function yields an average time constant of ~ 21 minutes 
and a x^/dof of 107.7/98, consistent with the hypoth- 
esis of a random arrival time of cosmic-ray muons, as 
expected. 



E. Expected Neutrino Signal and Background 
Rates 

Neutrino-induced muons from the H2O and rock sur- 
rounding the SNO detector were simulated in the manner 
outlined in the previous sections. A zenith angle cut of 
cos^zonithf: 0.4 was imposcd to reject cosmic-ray muons 
from our neutrino-induced signal. Under the assump- 
tions of the Bartol |22j atmospheric neutrino flux and no 
oscillations, SNO expects a total of 138.4 ± 7.3 neutrino- 
induced events per year passing all cuts. A full break- 
down of the expected signal contribution is shown in Ta- 
ble Iml 

The efficiency for reconstructing these signal events is 
not as high as that for primary cosmic-ray muons be- 
cause some of neutrino-induced events stop within the 
detector volume. The total efficiency is defined as the 
ratio between the number of through-going cosmic rays 
that reconstruct with an impact parameter less than 830 
cm that pass all cuts versus the number of through-going 
muon events with a generated impact parameter less than 
830 cm. If the muon is genuinely through-going (exits the 
fiducial area of the detector), the total efficiency is 98.0%, 
based on Monte Carlo studies. 

SNO also has a small acceptance for neutrino-induced 
muons whose interaction vertex resides inside a fidu- 
cial volume defined by the 830-cm radius. Most of 
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TABLE III: Summary of Monte Carlo expected signal and background rates contributing to the neutrino-induced muon analysis, 
after all cuts, for the full zenith angle range of —1 <cos0zonith< 0.4 and the u noscillated region of <cos 6'zcnith< 0.4. Errors 
include full systematic uncertainties assuming no correlations (see Table VII for more details). Neutrino induced interaction 



rates assume no oscillations. The last entry in the table shows the measured muon rate passing all cuts. 



Source 


Rate (yr"^) 


Zenith Range 


— 1 < COS6'zonith< 0.4 


< COS6'zenith< 0.4 


Through-going rock interactions 
Through-going water interactions 
Internal interactions 
Internal Ve interactions 
Total Signal 


iz4.4 ± 0.5 

9.0 ±0.5 

3.1 ±0.8 
1.9 ±0.3 

138.4 ±7.3 


,10 1 o o 

43.2 ± 2.3 
2.8 ±0.2 
1.0 ±0.3 
0.7±0.1 

47.8 ±2.5 


Cosmic ray /i 
Instrumental contamination 
Total Background 


1.1 ± 1.2 
0.3 ±0.2 
1.4 ± 1.2 


1.1 ± 1.2 
0.1 ±0.2 

1.2 ± 1.2 


Total Expected Rate 


139.8 ±7.4 


49.0 ±2.8 


Detected Rate 


152.7 


59.7 



these events are removed by the energy loss cut. From 
Monte Carlo studies, contamination of 3.1 ± 0.8 con- 
tained I'^-events per year is expected in the data. As 
these events also depend on the flux and neutrino oscil- 
lation parameters, they are included as part of the final 
signal extraction. A small number of internal neutrino 
events also come from v,, interactions which reconstruct 
as through-going muons. The rate of these events is 
1.9 ± 0.3 events/year. The cosmic-ray muon background 
passing all cuts is estimated to be 1.1 ± 1.2 events per 
year. Finally, a negligible amount of instrumental back- 
grounds are expected to contaminate the muon signal. 
The majority of such instrumentals are due to burst ac- 
tivity present in the detector. A bifurcated analysis com- 
paring the high-level cuts against the low-level cuts is per- 
formed so as to determine the amount of contamination 
of these instrumental events in our data jUJHS]. In addi- 
tion, events that are explicitly tagged as burst events are 
used to test the cut effectiveness in removing instrumen- 
tal contamination. Both tests predict an instrumental 
background contamination rate of 0.3 ± 0.2 events per 
year. 



V. FLUX AND OSCILLATION RESULTS 

A. Cosmic-Ray Muon Flux 

In order to minimize the possibility of introducing bi- 
ases, a two-tier blind analysis procedure is employed. 
First, only a fraction (~ 40%) of the data was open for 
analysis. Second, a fraction of muon events was removed 
from the data set using a zenith angle-dependent weight- 
ing function unknown to the analyzers. Only after all 
fitter and error analyses were completed were both blind- 
ness veils lifted. 



A total of 76749 muon candidates passing all se- 
lection cuts are reconstructed with a zenith angle of 
0.4 <cos6'zonith< 1 for the 1229.30-day dataset. The 
data collected corresponds to an exposure of 2.30 x 10^^ 
cm^ s. The total measured cosmic-ray muon fiux at SNO, 
after correcting for acceptance, is (3.31 ± 0.01 (stat.) ± 
0.09 (sys.)) X 10-1° /x/s/cm^, or 62.9 ± 0.2 muons/day 
passing through a 830-cm radius circular fiducial area. 

One can define the vertical muon intensity per solid 
angle by the expression: 

1 ^' 

cos (0,,,) - ^ ^ ^ ^ (3) 

where Ni is the number of events in a given solid angle 
bin f2.i and zenith angle O^^i, L is the livetime of the mea- 
surement, e is the detection efficiency for through- going 
muons, and A is the fiducial area. Given the fiat over- 
burden, it is possible to express Equation [3] in terms of 
the slant depth, xsno- To compare to other vertical fiux 
measurements, SNO rock can be corrected to standard 
rock, CaCOa, using the relation: 

xstd = 1.015XSNO ± ^^IF (4) 

4 X lO^m.w.e 

where xgNO is the slant depth expressed in meters water 
equivalent. There exists an additional ±1% model uncer- 
tainty in converting from SNO to standard rock which is 
estimated from differences that arise between the MU- 
SIC and PROPMU energy loss models. Flux values for 
slant depths ranging from 6 to 15 km water equivalent 
are presented in Table |IV[ 

The attenuation of the vertical muon intensity as a 
function of depth can be parameterized by: 
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TABLE IV: Intensity for standard (CaCOa) rock as a function of slant depth (in meters water equivalent) for muons passing 
all cuts and which reconstruct with cos ^zenith > 0.4. Only statistical errors are shown. 



Slant depth 


Events 


Intensity (standard rock) 


(meters w.e.) 




(cm~^ s^"' sr^""") 


6225 


4203 


(3.71 ±0.53) X 10"^° 


6275 


3905 


(3.47 ±0.50) X 10"^° 


6325 


3576 


(3.20 ±0.46) X 10"^° 


6375 


3371 


(3.05 ±0.44) X 10"^° 


6425 


3238 


(2.95 ±0.43) X 10"^° 


6475 


3000 


(2.75 ±0.40) X 10"" 


6525 


2737 


(2.53 ±0.37) X 10"^" 


6575 


2598 


(2.42 ±0.36) X 10"^" 


6625 


2369 


(2.23 ± 0.33) X 10"^° 


6675 


2182 


(2.07 ± 0.31) X 10~^° 


6725 


2038 


(1.94 ±0.29) X 10"^° 


6775 


1911 


(1.84 ±0.28) X 10"^° 


6825 


1831 


(1.77 ±0.27) X 10""' 


6875 


1668 


(1.63 ±0.25) X 10""' 


6925 


1552 


(1.52 ±0.24) X 10"^° 


6975 


1377 


(1.36 ±0.21) X 10"^° 


7025 


1359 


(1.35 ±0.21) X 10-1° 


7075 


1247 


(1.25 ±0.20) X 10-1° 


7125 


1163 


(1.18 ±0.19) X 10-1° 


7175 


1111 


(1.13 ±0.18) X 10-1° 


7225 


1043 


(1.07 ±0.17) X 10"i° 


7275 


910 


(9.40 ± 1.50) X 10"ii 


7325 


897 


(9.33 ± 1.50) X 10-11 


7375 


830 


(8.69 ± 1.40) X 10-11 


7425 


790 


(8.33 ± 1.40) X 10-11 


7475 


758 


(8.05 ± 1.30) X 10-11 


7525 


683 


(7.30 ± 1.20) X 10"ii 


7575 


713 


(7.67 ± 1.30) X 10"ii 


7700 


2241 


(6.13 ± 1.00) X 10"ii 


7900 


1791 


(5.03 ± 0.86) X 10-11 


8100 


1378 


(3.97 ±0.69) X 10-11 


8300 


1097 


(3.24 ±0.58) X 10-11 


8500 


859 


(2.60 ± 0.47) X 10-11 


8700 


670 


(2.08 ±0.39) X 10"ii 


8900 


504 


(1.60 ±0.31) X 10-11 


9100 


444 


(1.44 ±0.28) X 10-11 


9300 


328 


(1.09 ±0.22) X 10-11 


9500 


257 


(8.7 ± 1.8) X IQ-i^ 


9700 


205 


(7.1 ± 1.5) X 10-1^ 


9900 


183 


(6.5 ± 1.4) X IQ-i^ 


10250 


291 


(4.3 ±0.9) X IQ-i^ 


10750 


166 


(2.6 ±0.6) X IQ-i^ 


11250 


100 


(1.6 ±0.4) X IQ-i^ 


11750 


61 


(1.0 ±0.3) X 10"" 


12250 


34 


(6.0 ± 1.8) X 10"i'^ 


12750 


31 


(5.7 ± 1.7) X 10-" 


13250 


14 


(2.7 ± 1.0) X 10-" 


13750 


10 


(2.0 ±0.8) X 10-" 


14250 


11 


(2.3 ±0.9) X 10-" 


14750 


7 


(1.5 ±0.7) X 10-" 


15250 


13 


(2.9 ± 1.1) X 10-" 
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/;(x.td) = /o(— )" e- 



(5) 



where /q is an overall normalization constant, and Xq 
represents an effective attenuation length for high-energy 
muons. The remaining free parameter, a, is strongly cor- 
related with the spectral index 7 in Equation |2] Results 
from fits of the vertical muon intensity as a function of 
depth for various values of these parameters are shown 
in Table |V] We perform fits whereby the parameter a 
is either fixed to what one would expect from the sur- 
face (a = 7 — 1 = 1.77) or allowed to float freely. The 
cosmic-ray data tends to prefer larger values of a than 
the expected value of 1.77. A comparison of SNO's muon 
flux to that measured in the LVD [46j and MACRO [47j 
is shown in Figure [9j In general, there exists tension be- 
tween the different data sets. Fits have been performed 
both with and without allowing the slant depth uncer- 
tainty to float within its uncertainty. The fits in both 
cases are nearly identical, with minimal change (< Icr) 
to the slant depth. The fits presented in Table fV] are with 
the slant depth constrained. 

To avoid some of the strong correlations between the 
three parameters listed in Equation |5] we also perform 
the fit using the following paramctrization: 



where e"" represents the muon fiux at the surface, oi is 



inversely proportional to the muon attenuation length, 
and 02 represents the deviation from the simple expo- 
nential model. Results from fitting to Eq. [6]are shown in 
Table IVll 



The systematic uncertainties of this measurement are 



summarized in Table VII Certain systematic errors for 
the cosmic-ray muon flux are in common with those 
for the neutrino-induced muon flux results, including 
livetime, impact parameter bias, and angular resolu- 
tion. Others are unique to the cosmic-ray muon flux. 
These include uncertainties in the rock density, the sur- 
face variation, the rock conversion model, muon strag- 
gling, instrumental backgrounds, and backgrounds from 
neutrino-induced events and multiple muons. This last 
background is estimated from events measured from the 
MACRO experiment |3H] . As the reconstruction for mul- 
tiple muon events in the detector is not well known, we as- 
sign a ±100% uncertainty on this potential background. 
These systematic uncertainties are included as part of the 
total error presented in Table [V] 



B. Atmospheric Neutrino Results 

We assume a model for the atmospheric neutrino flux, 
and flt for a total flux scaling factor as well as the atmo- 
spheric neutrino oscillation parameters. In these fits we 
use a two-neutrino mixing model: 



a>(L/i?., Am'), = $0 • [1 - sin^ 29 ■ sin^ ( ^^^^|^)] 



(7) 



where 9 is the neutrino mixing angle. Am? is the square 
mass difference in eV^, L is the distance traveled by the 
neutrino in km, is the neutrino energy in GeV, and 
<i>o is the overall normalization of the neutrino-induced 
flux. 



Although the signal uncertainty is dominated by statis- 
tics, systematic errors do have an impact on both the 
acceptance and zenith angle distribution of events. To 
account for distortions in the zenith angle spectrum, we 
generalize the x'^-puU technique (see |49] and references 
therein) to the case of a maximum likelihood analysis. 
This allows us to account for the smallness of statistics 
while still incorporating any correlations that may ex- 
ist between different systematic error contributions. An 
extended likelihood function is constructed using the fol- 
lowing equation: 



Jvbins Tydata 

Aotai = 2( ^ In (^) - (A^f ^ - A^f *^)); (8) 

i ^ 

where jSf'^'^^'^i^C) represents the number of data (Monte 
Carlo) events found in a given zenith bin i. To account for 
the effect of systematic errors on our likelihood contours, 
we perform a linear expansion of N^'~' with respect to 
a nuisance parameter a for each systematic uncertainty 
such that: 



Nr^N^P+Y^{^)Aa,^N^P{l+(3,Ad) (9) 

Note that we have used vector notation to denote a 
summation over all nuisance parameters. By expanding 
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TABLE V: Results from the SNO fit to the vertical muon intensity for cos &zcnith> 0.4 using Equation[5] The fits were performed 
either using only SNO data with the a parameter allowed to float, with the a parameter fixed to the value predicted from the 
surface flux of Eq. [2] (a = 7 - 1 = 1.77), or combined with LVD [46] and MACRO t47j cosmic ray data. Symbols in the table 
are as defined in the text. The errors reported are a combination of statistical and systematic uncertainties on the fiux and 
slant depth. 



Dataset 


/o 

(10"'' cm-2 s"^ sr"^) 


(km w.e.) 


Q 


xVdof 


SNO only 
SNO only 
SNO + LVD + MACRO 


1.20 ±0.69 
2.31 ±0.32 
2.16 ±0.03 


2.32 ±0.27 
1.09 ±0.01 
1.14 ± 0.02 


5.47 ±0.38 

1.77 
1.87 ±0.06 


34.2 / 44 
111.0 / 45 
230.2/134 



TABLE VL Results from the SNO fit to the vertical muon intensity for cosSzenith> 0.4 using Equation |6] Fits shown using 
only SNO data, or combined with LVD [46] and MACRO [47] cosmic ray data. Symbols in table are as defined in the text. 
The errors reported are a combination of statistical and systematic uncertainties on the flux and slant depth. 



Dataset 


(cm~'^ s^-"" sr^"') 


ffli 

(m.w.e.)"^ 


ai 

(m.w.e.)-^ 


xVdof 


SNO only 
SNO ± LVD ± MACRO 


(4.551°;?°) X 10-" 
(1.97 ±0.06) X lO"'' 


(-1.75 ±0.06) X 10-3 
(-1.55 ±0.01) X 10-3 


(3.9 ±0.3) X 10^* 
(2.78 ±0.08) X 10"** 


41.6 / 44 
230.8 / 134 




1 n-14l I \ I \ \ I \ \ \ \ I \ \ I \ I I \ I I I I I 

4000 6000 8000 10000 12000 14000 

Slant Depth (km w.e.) 

FIG. 9: The flux of cosmic-ray muons that pass all cuts as 
a function standard rock depth. SNO data (filled circles) 
shown with best global fit intensity distribution (dashed line) 
and data from LVD [46] (empty circles) and MACRO [47] 
(triangles) detectors using Eq. [5] Global fit range extends to 
13.5 kilometers water equivalent, beyond which atmospheric 
neutrino-induced muons start to become a significant fraction 
of the signal. 



the logarithmic term to second order and minimizing the 
likelihood function with respect to each nuisance param- 
eter, one finds an analytical expression [50j : 
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-A/'bins Tydata 

Aotal = 2( ^ In (-^) - «C _ ^data)) _ Ad^.^S^ Aa^,^ (10) 



where AQ?^j„ represents the minimized nuisance param- 
eter: 



iVbins 

Aa„„ = ( {Nt^"^ - KfmS-^ (11) 

i 

and the matrix 5^ is defined as: 



ing parameters 6 and Am^ in Equation [t] A max- 
imum likelihood fit is performed to find the best fit 
points, as outlined above. If all parameters are al- 
lowed to float, one finds a fiux normalization value of 
<i>o =1.22 ±0.10 and best fit neutrino oscillation parame- 
ters of Ato^ of 2.6 X lO"'^ eV^and maximal mixing. These 



(12) 



Here, ct^^ is the diagonal error matrix whose entries 
represent the size of the systematic error constraints. 
As long as the contribution from the systematic errors 
is small, the above formalism provides a very efficient 
method for evaluating the effect of systematic errors 
while also incorporating constraints from the data. A 
total of six systematic uncertainties are fit using this 
method; five of which (axial mass, quasi-elastic cross- 
section, resonance cross-section, deep inelastic scatter- 
ing, and energy loss modeling) have explicit zenith angle 
dependencies, while the last is fiat with respect to the 
zenith distribution. This uncertainty is a combination of 
all of the remaining systematic errors and is fit as an over- 
all normalization error. A summary of all the systematic 



errors is shown in Table WJu 

Figure [lO] shows the zenith angle distribution for 
neutrino-induced muons. A total of 514 events are 
recorded with —1 <cos 6'zonith< 0.4 in the 1229.30 days 
of livetime in this analysis. For neutrino-induced 
events near the horizon (cos^zenith between and 0.4), 
201 events are observed. Given the current measurements 
of the atmospheric oscillation parameters, the neutrino- 
induced flux is unaffected by oscillations in this latter 
region and therefore is a direct measurement of the atmo- 
spheric neutrino flux, particularly at high energies. The 
corresponding neutrino-induced through-going muon flux 
below the horizon (cos0zenith< 0) and above the horizon 
(0 <cos 6lzenith< 0.4) are 2.10 ± 0.12(stat.) ± 0.08(sys.) x 
10-i3cni-2s-isr-i and 3.31 ± 0.23(stat.) ± 0.13(sys.) x 
10~^^cm~^s^^sr~^, respectively. 

From the measured zenith angle distribution, we can 
extract the flux normalization $0 and the neutrino mix- 
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FIG. 10: The distribution of through-going neutrino-induced 
muons that pass all cuts as a function of zenith angle. 
Data (crosses) are shown with the best-fit MC spectra of 
($o,sin^ 261, Am^) = (1.22 ±0.10, 1.00, 2.6 x 10"^ eV^) (solid 
box) and prediction with no neutrino oscillation and a best 
fit normalization of $0 = 1.09±0.08 (hashed box). The back- 
ground due to cosmic-ray muons is shown in the dashed line. 
The zenith angle cut is indicated in the figure. 



results are with respect to the Bartol three-dimensional 
atmospheric flux model and the cross-section model im- 
plemented in NUANCE described in Section|l!!][22 . The 
zenith angle spectrum is consistent with previously mea- 
sured neutrino oscillation parameters. One can also look 
at SNO's sensitivity on the atmospheric flux <I>o by in- 
cluding existing constraints on the atmospheric neutrino 
oscillation parameters from the Super-Kamiokande |TT] 
(Am2,sin2 26lsK) = (2.i;o.4 x 10"^ eV^ 1.000 ± 0.032) 
and MINOS ill [SI] (Am^jNos = (2-43 ± 0.13) x 
10~^ eV^) neutrino experiments. The likelihood function 
in Eq. [TOlis altered to the following: 



r _r , , Am^-Am|K ^2 , / ^"i^ - ^"^minos ^2 , , sin^ 20 - sin^ 26Isk ,2 (.0. 

•i-constrainod — i-total + ) +\ ) + ( ) ('-•j) 



CTAm^.SK CTAm^.MINOS Cg.SK 
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TABLE VII: Summary of systematic errors for the neutrino-induced and cosmic-ray muon flux measurements. A dagger (f) 
indicates that the systematic error only affects the cosmic-ray muon intensity fit to Eq. [5] and is not included in the total 
systematic error summation below. The total error in the table is determined from the fit including correlations and does not 
equal to the quadrature sum of the individual components. 



Systematic Error 


Variation 


!^;j-induced muon flux error 


Cosmic-ray muon flux error 


Detector 


Detector Propagation Model 


Various 


±0.3% 


±0.3% 


Angular Resolution 


±0.6° 


±0.1% 


±0.1% 


Energy Loss Model 


±5% 


±2.5% 


±0.2% 


Impact Bias/Shift 


±4.0 cm 


±1.2% 


±1.0% 


Impact Resolution 


±8.5 cm 


±0.07% 


±0.07% 


Livetime Clock 


±2600 s 


±0.002% 


±0.002% 


PMT Charge Model 


±10% 


±1.0% 


±0.05% 


Fisher Discriminant Cut 


±5% 


±2.1% 


±0.37% 


Total Detector Model 




±3.7% 


±1.1% 


Neutrino Cross-Section Model 


Axial Mass 


±0.15 GeV 


±1.1% 


N/A 


Quasi-Elastic 


±10% 


±0.8% 


N/A 


Resonance 


±20% 


±1.9% 


N/A 


Deep Inelastic 


±3% 


±2.1% 


N/A 


Total Cross-Section Model 




±3.1% 


N/A 


Muon Propagation Model 


Rock Density(f ) 


±0.05 g/cm^ 


±0.3% 


(t) 


Conversion Model(f) 


±1% 


N/A 


(t) 


Surface Variation(f ) 


±50 m 


N/A 


(t) 


Transport Model 




±2% 


N/A 


Time/Seasonal Variation 




±1% 


±2.2% 


Total Propagation Model 




±2.2% 


±2.2% 


Backgrounds 


Instrumental 


0.3 ± 0.2 events yr"^ 


±0.2% 


< 0.1% 


Cosmic ray fi 


0.6 ± 1.1 events yr~^ 


±0.8% 


N/A 


i^^j-Induced 


45.8 ± 2.3 events yr"^ 


N/A 


±0.2% 


Multiple Muons 


±100% 


« 1% 


±1% 


Total Background Error 




±0.8% 


±1% 


Total Systematic Error 


±4.8% 


±2.7% 


Statistical Error 


±8.5% 


±0.4% 



The constraint reduces the uncertainty on the overall 
atmospheric neutrino flux normalization to 1.22 ± 0.09. 
The 68%, 95% and 99.73% confidence level regions for 
the parameters as determined by the fits are shown in 
Figure 11. The scenario of no neutrino oscillations by 
using SNO-only data is excluded at the 99.8% confidence 
level. 



VI. SUMMARY 

The Sudbury Neutrino Observatory experiment has 
measured the through-going muon flux at a depth of 



5890 meters water equivalent. We find the total muon 
cosmic-ray fiux at this depth to be (3.31 ± 0.01 (stat.) ± 
0.09 (sys.)) X 10^^° ^/s/cTT?. We measure the through- 
going muon flux induced by atmospheric neutrinos. The 
zenith angle distribution of events rules out the case of 
no neutrino oscillations at the 3cr level. We measure the 
overall flux normalization to be 1.22±0.09, which is larger 
than predicted from the Bartol atmospheric neutrino flux 
model but consistent within the uncertainties expected 
from neutrino flux models. This is the first measure- 
ment of the neutrino-induced flux above the horizon in 
the angular regime where neutrino oscillations are not an 
important effect. The data reported in this paper can be 
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FIG. 11: The 68% (blue), 95% (red), and 99.73% (black) confidence level contours for the i/^ atmospheric neutrino oscillation 
parameters based on the muon zenith angle distribution for cos0zcnith< 0.4. The plots show the SNO-only contours for fiux 
normalization versus mass splitting (top left), SNO-only mass splitting versus mixing angle (top right), SNO-only contours for 
flux normalization versus mixing angle (bottom left) and the flux normalization versus mixing angle including constraints from 
the Super-K and MINOS neutrino oscillation experiments (bottom right) [111 I13j . 



used to help constrain such models in the future. 
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